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Abstract 

We explore the combined physics potential of T2K and NOz/A in light of the moderately large 
measured value of #13. For sin 2 26*13 = 0.1, which is close to the best fit value, a 90 % C.L. evidence 
for the hierarchy can be obtained only for the combinations (Normal hierarchy, —170° < 6c p < 0°) 
and (Inverted hierarchy, 0° < 5cp < 170°), with the currently planned runs of NOi'A and T2K. 
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are increased by 50% and those of T2K are doubled. Such an increase will also give an allowed 
region of 5cp around the its true value, except for the CP conserving cases 5cp = or ± 180°. We 
demonstrate that any measurement of 5c p is not possible without first determining the hierarchy. 
We find that comparable data from a shorter baseline (L ~ 130 km) experiment will not lead to 
any significant improvement. 
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I. INTRODUCTION 



Neutrino physics has entered a phase of precision measurements. During the past few 
years, the following precise measurements of neutrino parameters have been made with high 
intensity sources: 

• The smaller mass-squared difference A21 = m\ — m\ is measured by KamLAND [1] 
while the precision on # 12 is controlled by the solar experiments [2] . Global analysis of 
all the data, in the three flavour oscillation framework, gives A 21 = (7.6 ± 0.2) x 1CT 5 
eV 2 and sin 2 12 = 0.312 ± 0.016 [3]. 

• MINOS [4] experiment has measured the magnitude of the mass-squared difference in 
the survival probability. The precision on # 23 is controlled by atmospheric neutrino 
data [5]. Global analysis gives two distinct values of A 31 depending on whether it is 
positive [which is the case for normal hierarchy (NH)] or negative [which is the case 
for inverted hierarchy (IH)]. The ranges are A 31 (iVif) = (2.45 ± 0.09) x 10 -3 eV 2 and 
A 31 (IH) = (-2.31 ± 0.09) x lO" 3 eV 2 with sin 2 6 23 = 0.51 ± 0.06 for both cases [3]. 

• The global fits to data from the accelerator experiments T2K [6] and MINOS [7] and 
the reactor experiments DChooz [8], Daya Bay [9] and RENO [10] have determined 
#13 to be non-zero at 5a level, with the best fit very close to sin 2 2# 13 ~ 0.1 [11, 12]. 

We expect the following improvements in precision during the next few years. 

• Very high statistics data from T2K [13] and MINOS [4] experiments will improve the 
precision on |A 31 | and sin 2 2# 23 to a few percent level. 

• Reactor experiments are taking further data [14-17]. The survival probability at 
these reactor experiments is sensitive only to the mixing angle #i 3 and hence they can 
measure this angle unambiguously. By the time they finish running (around 2016), we 
estimate that they should be able to measure sin 2 20i 3 to a precision of about 0.005. 

In light of these current and expected near future measurements, the next goals of neutrino 
oscillation experiments are the determination of neutrino mass hierarchy, detection of CP 
violation in the leptonic sector and measurement of Sep- These goals can be achieved by high 
statistics accelerator experiments measuring — >■ v e and — » u e oscillation probabilities. 
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Among such experiments, T2K is presently taking data and NOz/A is under construction 
and is expected to start taking data around 2014. All other experiments, capable of making 
these measurements, are far off in future. In this paper, we study the combined ability of 
T2K and NOz/A to achieve the above goals. 

In the above discussion, we have two different magnitudes for A31 for the two hierarchies 
because the mass-squared difference measured in survival probability is not A31 but is an 
effective one defined by [18, 19] 

Am^ = A 3 i - (cos 2 9 12 - cos 5 C p sin 9 13 sin 29i 2 tan #23) A 2 i- (1) 

Accelerator experiments, such as MINOS and T2K, measure the magnitude of the above 
quantity. But the magnitudes of A31 will turn out to be different for A31 positive (NH) and 
A31 negative (IH). 

II. SIMULATION DETAILS 

Before discussing various physics issues, we discuss the details of our simulation. We do 
this because we will illustrate various points through the means of simulation. 

We use the software GLoBES [20, 21] for simulating the data of T2K, NOz/A and an 
envisaged short baseline experiment from CERN to Frejus (C2F), which is a scaled down 
version of MEMPHYS [13, 22-32]. Various details of these experiments and their charac- 
teristics, especially the signal and background acceptances, are given in Table I. The basic 
properties of NOz/A are taken from Ref. [24] and of T2K are taken from Ref. [13]. The 
efficiencies for each of the experiments are taken from GLoBES [20, 21]. The background 
errors consist of errors in flux normalization (norm) and in spectrum (tilt). 

We have kept the solar parameters A21 and 9 12 fixed at their best fit values throughout 
the calculation. We have taken the central values of |A 31 | and #23 to be their best fit values. 
We took a (sin 2 2^23) = 0.02 and a (| A 31 |) = 0.03 x (| A 31 |), because of the precision expected 
from T2K. We have done computations for various different values of sin 2 2#i3 in the range 
0.05 — 0.2 [11, 12]. We took o (sin 2 26 l i3) = 0.005 which is the final precision we can hope 
for from the reactor experiments. The value of the CP-violating phase 8c p is varied over its 
entire range -180° to 180°. 
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Characteristic 


NOz^A 


T2K 


C2F (assumed) 


Baseline 


812 km 


295 km 


130 km 


Location 


Fermilab - Ash River 


J-PARC - Kamioka 


CERN - Frejus 


Beam 


NuMI beam 0.8° off - 


JHF beam 2.5° off - 


SPL superbeam 




axis 


axis 




Beam power 


0.7 MW 


0.75 MW 


0.75 MW 


Flux peaks at 


2 GeV 


0.6 GeV 


0.35 GeV 


P^e 1st Osc. Maximum 


1.5 GeV 


0.55 GeV 


0.25 GeV 


Detector 


TASD, 15 kton 


Water Cerenkov, 22.5 


Water Cerenkov, 22.5 






kton 


kton 


Runtime (years) 


3 in v + 3 in D 


5 in v 


3 in v + 3 in v 


Signal 1 (acceptance) 

O \ ST / 


v e appearance(26%) 

c ST ST \ v / 


u e appearance(87%) 


v e appearance(71%) 

c sr sr \ " / 


Signal 1 error 


5%, 2.5% 


2%, 1% 


2%, 0.01% 


(norm., tilt) 








Background 1 (acceptance) 


mis - id muons/anti 


mis - id muons/anti 


mis - id muons/anti 




- muons(0.13%), NC 


- muons, NC events, 


- muons(0.054%), NC 




events(0.28%), Beam 


Beam f e /P e (binned 


events(0.065%), Beam 




i/ e /P e (16%) 


events from GLoBES 


z/ e /P e (70%) 






[20, 21]) 




Background 1 error 


10%, 2.5% 


20%, 5% 


2%, 0.01% 


(norm., tilt) 








Signal 2 (acceptance) 


v e appearance(41%) 

c sr ST \ v / 


u e appearance(87%) 

c sr sr V / 


v e appearance(68%) 

c sr sr \ v j 


Signal 2 error 


5%, 2.5% 


2%, 1% 


2%, 0.01% 


(norm., tilt) 








Background 2 (acceptance) 


mis - id muons/anti 


mis - id muons/anti 


mis - id muons/anti 




- muons(0.13%), NC 


- muons, NC events, 


- muons(0.054%), NC 




events(0.88%), Beam 


Beam v e lv e (binned 


events(0.25%), Beam 




j/ e /P e (33.6%) 


events from GLoBES 


^e/Pe(70%) 






[20, 21]) 




Background 2 error 


10%, 2.5% 


20%, 5% 


2%, 0.01% 


(norm., tilt) 
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TABLE I: Properties of various long baseline experiments considered in this paper. 



We compute statistical 



Xlt as 



true 




(2) 



where Nj^ rue is the event distribution for true hierarchy and some fixed true value of Sep- 
jytest j g ^ e event distribution with the test hierarchy either true or wrong and a varying 
test value of 6c p as inputs. The index i runs over the number of energy bins. The final \ 2 
is computed including the systematic errors, described in Table I, and the priors on IA31I, 
sin 2 2^23 and sin 2 26 l i3. The prior on sin 2 2^13 effectively takes into account the data due to 
reactor neutrino experiments. 

In the following we consider two kinds of plots both of which are shown as contours in 
the sin 2 29 13 -5cp plane. 

• Hierarchy exclusion plots: These are plotted in the plane of true values of sin 2 26\ 3 -8cp- 
The contours in these plots define the line \ 2 = 2.71. In computing this x 2 > we have 
marginalized over the parameter ranges described above. For all sets of parameter 
values to the right of the contour, the wrong hierarchy can be ruled out at 90% C.L. 

• Allowed region plots: These are plotted in the plane of test values of sin 2 29i 3 -Scp- 
The contours in these plots are defined by x 2 = 4.61. The region enclosed by them 
is the set of allowed values of sin 2 29is-5cp at 90% C.L. for a given set of neutrino 
parameters. 

Throughout this paper, the phrase "hierarchy determination" implies 90% C.L. evidence 
for hierarchy. 

III. HIERARCHY DETERMINATION WITH P Me 

The — > v e channel is sensitive to a number of neutrino parameters and hence is the 
most sought after in the study of neutrino oscillation physics using long baseline experiments. 
In the presence of matter, the — > v e oscillation probability, expanded perturbatively in 
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the small mass-squared difference, A 2 i is given by [33-35] 

9 „ 9 „ sin 2 A(l - A) 
P -> v e ) = P^ = sin 2 2fl 13 sm 2 23 



+a cos 6\3 sin 2#i 2 sin 26*i3 sin 2# 23 cos( A + Sep 



1 - A) 2 

sin A^sin A(l - A) 



A I- A 



+a 2 sin 2 26 12 cos 2 0^ cos 2 6 23 Sm (3) 

A 2 

where A = A 31 L/4P, A = A/A 31 , a = A 2 i/A 31 . A is the Wolfenstein matter term [36] and 
is given by A(eV 2 ) = 0.76 x l(T 4 p (gm/cc)P(GeV). 

For NH A 3 i is positive and for IH A 3 i is negative. The matter term A is positive for 
neutrinos and is negative for anti-neutrinos. Hence, in neutrino oscillation probability, A is 
positive for NH and is negative for IH. For anti-neutrinos, A is negative for NH and positive 
for IH and the sign of 8c p is reversed. The presence of the term A in P^ e and in P^ g makes 
them sensitive to hierarchy. The longer the baseline of an experiment, the greater is the 
sensitivity to hierarchy because, P^ e peaks at a higher energy for longer baseline and the 
matter term is larger for higher energies. 

As can be seen from Eq. (3), P Me is dependent on #i 3 , hierarchy amd 8c p in addition to 
other well determined parameters. A measurement of this quantity will not give us a unique 
solution of neutrino parameters but instead will lead to a number of degenerate solutions 
[37-40]. Since 13 is measured unambiguously and precisely [8-10], degeneracies involving 
this parameter are no longer relevant. Only hierarchy-^cp degeneracy has to be considered. 
This degeneracy prevents any one experiment from determining hierarchy and 8c p, leading 
to the need for data from two or more long baseline experiments [22, 41-44]. 



A. Hierarchy- 5c p degeneracy for NOz^A 

First we consider the hierarchy determination capacity of NOz/A alone because the matter 
term and the hierarchy dependence is the largest for this experiment, due to the flux peaking 
at higher energy. In Fig. 1 (left panel), we have plotted P Me vs E for both NH and for IH 
for NOz/A baseline of 812 km. The bands correspond to the variation of 8cp from —180° 
to +180°. The values of P^ e are, in general, higher for NH and lower for IH. This is a 
straightforward consequence of the A dependence of P^ e . Further, we note that for both 
NH and IH, the value of 8c p = +90° gives the lowest curve in the band and the value 
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of Sep = —90° gives the highest curve in the band. This behaviour can also be easily 
understood from Eq. (3). At the oscillation maximum, A ~ 90°. Hence cos(A + Sep) is 
+1 for Sep = —90° and is —1 for Sep = +90°. As can be seen from the figure, there is an 
overlap of the bands for (NH, Sop — +90°) and (IH, Sep — —90°). Hence, if the measured 
probability comes to be these values, then we have two degenerate solutions. In Fig. 1 (right 
panel), we have plotted the corresponding anti-neutrino probabilities. P^ g is higher for IH 
and lower for NH as a consequence of the reversal of A sign. Since Sep sign is reversed for 
anti-neutrinos, here Sep = +90° defines the upper curves and Sep = —90° defines the lower 
curves. Here again there is an overlap between (NH, Sep — +90°) and (IH, Sep — —90°) so 
we get the same degenerate solutions as the neutrino case. 
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FIG. 1: (colour online) P^e (left panel) and P^ g (right panel) bands for NOz/A for sin 2 2(9i3 = 0.1 



From Fig. 1, we can define the concept of favourable half plane for each hierarchy. Suppose 
NH is the true hierarchy. If Sep is in the lower half plane (—180° < Sep < 0°, LHP) then 
all the curves for P^NH^Scp) lie much above the set of curves for P^ e (IH,Scp)- In the 
case of anti-neutrinos, P^NH, Sep) will be much below P^IPi, Sep)- In such a situation, 
the data from NOz/A alone can determine the hierarchy. Therefore we call the LHP to be 
the favourable half-plane for NH. Similar arguments hold if IH is true hierarchy and Sep is 
in the upper half plane (UHP). So UHP is the favourable half plane for IH. Thus, if nature 
chooses one of the following two combinations (NH, LHP) or (IH, UHP), then NOz/A , by 
itself, can determine the hierarchy. 
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The separation between the set of curves P^ e (NH,Scp) and P^IH^Scp) also depends 
on 8 13 . The two sets have more overlap for smaller values of (9 13 but become more separated 
for larger values of #13. This is illustrated in Fig. 2, showing P^ e vs E, for a lower and higher 
value of sin 2 2#i3 . It is easier to determine the hierarchy if the separation between the curves 
is larger, that is if #13 is larger. 
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FIG. 2: (colour online) P^ e bands for NO^A for sin 2 26*13 = 0.05 (left panel) and 0.15 (right panel) 



The favourable and unfavourable half planes for a particular hierarchy can also be defined 
from Eq. (3), where the Sep dependence occurs purely in the form cos(A + Sep)- If NH is 
the true hierarchy, A 90° around the probability maximum. Then, the Sep dependent 
term increases P^ if Sep is in the LHP and decreases it if Sep is in the UHP. Hence a 
cleaner separation from P^IH^Scp) can be obtained only if Sep is in the LHP. If IH is 
the true hierarchy, A ps —90°. Then P^ e is reduced, and moved away from P^NH, Sep) if 
Sep is in the UHP. Thus UHP forms the favourable half plane for IH, whereas LHP is the 
favourable half plane for NH. Even if we use the anti-neutrino oscillation probabilities, the 
same considerations will hold. Therefore, the same relation between hierarchy and half-plane 
holds for both neutrino and anti-neutrino data. 

We plot the hierarchy discrimination ability of NOi^A in Fig. 3. We see that, for 
sin 2 2^13 = 0.1, the hierarchy can be determined at 90 % CL. for the following two com- 
binations: (NH, -170° < S C p < -10°) or (IH, 10° < S C p < 170°). The statistics for the 
experiment are not quite enough to determine the hierarchy for the whole favourable half 
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plane for this value of #13. If sin 2 29\s = 0.12, then the hierarchy can be determined for the 
whole favoured half plane. It was shown in Ref. [45] that NOz/A can determine the hierarchy 
for 45 % of the 5cp range for sin 2 29 13 = 0.1. 
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FIG. 3: (colour online) Hierarchy exclusion plots for NO^A for Zv+W running when NH is true (left panel) 
and when IH is true (right panel) 



For smaller values of sin 2 2#i3, one needs larger statistics to determine the hierarchy for 
the whole favourable half plane. This is illustrated in Fig. 4. With 1.5 times the presently 
projected statistics of NOi/A , one can determine the hierarchy for the whole of the respective 
favourable half planes, for both NH and IH, for sin 2 2^13 = 0.1. Similar conclusions were 
obtained earlier in Ref. [46]. If 5cp happens to be in the unfavourable half plane, even 
tripling of statistics leads to hierarchy determination only for a very small range of 5c p- 
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True:NH / Test:IH, 90% C.L. True:IH / Test:NH, 90% C.L. 




FIG. 4: (colour online) Hierarchy exclusion plots for NOi^A with boosted statistics for 3v+3iy running 
when NH is true (left panel) and when IH is true (right panel). 

B. Resolving the hierarchy-^cp degeneracy with T2K 

As we demonstrated in the previous subsection, NOi^A alone can't determine the hier- 
archy if nature chooses one of the unfavourable combinations (NH, UHP) or (IH, LHP). In 
this subsection, we explore how data from T2K can help in resolving this problem. Since 
the baseline of T2K is smaller, the probability peaks at a lower energy and hence the flux is 
designed to peak at a lower energy. Therefore the matter term A is much smaller for T2K. 
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True:NH / Test:IH, 90% C.L. 



True:IH / Test:NH, 90% C.L. 
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FIG. 5: (colour online) Hierarchy exclusion plots for NOz^A + T2K with nominal statistics when NH is 

true (left panel) and when IH is true (right panel). 



In Fig. 5, we plot the combined hierarchy exclusion capability of NOi^A and T2K. From 
this figure we see that, for sin 2 26* 13 < 0.1, hierarchy determination is not possible for any 
5cp in the unfavourable half-plane, Hence, in our example, we assume that the statistics of 
NOi^A are 50 % more than the nominal value and those of T2K are twice the nominal value. 

We illustrate the effect of T2K data on hierarchy determination by a set of examples. 
First we assume that NH is the true hierarchy and the true value of Sep = 90°, i.e. in the 
unfavourable half plane. In such a situation, NOz/A data gives two degenerate solutions in 
the form of (NH, 5cp ~ 90°) and (IH, 5qp in LHP), as shown in Fig. 6 (left panel). 
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1.5*N0vA, 90%C.L. 



1.5*NOvA+2*T2K, 90% C.L. 
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FIG. 6: (colour online) Allowed sin 2 26 13 -6 C p plots for 1.5*NOA (left panel) and 1.5*NOA + 2*T2K 
(right panel) with true sin 2 29i 3 = 0.1 and true Sep = 90°. 
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1.5*NOvA+2*T2K, 90% C.L. 
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FIG. 7: (colour online) Allowed sin 2 29 13 -5 C p plots for 1.5*NOi/A (left panel) and 1.5*NOi/A + 2*T2K 

(right panel) with true sin 2 2#i 3 = 0.1 and true Sep = —90°. 

But, the addition of T2K data almost rules out the (IH,LHP) solution, seen in the right 
panel of Fig. 6. It is true that a very small part of the allowed region is left behind. But, 
comparing the two panels of Fig. 6, we see that the addition of T2K data reduces the 
allowed NH region only by a small amount whereas the allowed IH region is drastically 
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reduced. This gives a very strong indication of which hierarchy is correct. Thus the data 
of NOz/A in conjunction with that of T2K can effectively discriminate against the wrong 
hierarchy. This holds true for the case of IH being the true hierarchy with 5c p in LHP, 
illustrated in Fig. 7. Figs. 6 and 7 are similar to figures 2 and 3 of Ref. [45], which are done 
for the same 5c p values. Those figures also show the large shrinkage of the wrong hierarchy 
solution, with the addition of T2K data. In the following, we will demonstrate that this 
feature occurs for all values of 5<7p(true) in the unfavourable half-plane. 

A theoretical analysis of the hierarchy-^p degeneracy resolution, with data from NOz/A 
and T2K, was done in Ref. [47]. To keep the arguments simple, first it was assumed that 
023 is maximal and that sin 2 2#i3 is measured accurately by the reactor experiments. In such 
a situation, given a probability measurement, there exist two degenerate solutions: (correct 
hierarchy, correct 5cp) and (wrong hierarchy, wrong Sep)- In Ref. [47], it was shown that, 
for a given experiment, [sin(correct 5cp) —sin (wrong Sep)} is proportional to the matter term 
A for that experiment. For T2K, this difference is small and is about 0.7 for sin 2 2# 13 = 0.1. 
For NOz/A it is three times larger. Therefore, the wrong 5c p values for T2K data and for 
NOz/A data are widely different. A combined analysis of data from T2K and NOz/A will 
pick out the correct hierarchy and a range of 5c p around the correct value, provided the 
statistics from each experiment are large enough 

The above idea is illustrated below in Figs. 8 and 9. In Fig. 8, we have plotted x 2 vs 
<5cp(test) for various true values of Sep for NOz/A experiment. In the left panel, the true 
values of 5c p are all in LHP which is the favourable half-plane for NH. We find that, except 
for the CP conserving case of 5cp = —180°, all the \ 2 are above 9. Hence the wrong 
hierarchy can be excluded for most of the values of 5c p in the favourable half-plane. In 
the right panel, the true values of 5cp are all in UHP, which is the unfavourable half-plane 
for NH. And we find that in all cases, the \ 2 becomes nearly zero (except for 5cp = 0) 
for —120° < 5cp (test) < —60°. These are the degenerate (wrong hierarchy, wrong Sep) 
solutions mentioned above. Hence it is impossible to rule out the wrong hierarchy if true 
5cp is in the unfavourable half-plane. 
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1.5*N0vA True sin 26 13 : 0.1 



1 .5 *NOv A True sin 26 , 3 : 0.1 
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FIG. 8: (colour online) \ 2 vs - test Sep for 1.5*NOf A . Here true and test sin 2 26i 3 = 0.1. NH is true and 
IH is test. Different curves correspond to various true 8c p in lower half plane (left panel) and upper half 
plane (right panel). 



In Fig. 9, we have plotted x 2 vs 5cp{test) for various true values of Sep for T2K exper- 
iment. Once again, the left panel contains plots for 5c p in LHP and the right panel the 
plots for 5c p in UHP. From the left panel, we see that T2K can't rule the wrong hierarchy. 
This is to be contrasted with NOi^A case, where the wrong hierarchy is ruled out by NOi/A 
alone, if Sep is in the favourable half-plane. But, as we see below, T2K data is crucial for 
hierarchy discrimination, if 5c p is in the unfavourable half-plane. 
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FIG. 9: (colour online) \ 2 vs. test S CP for 2*T2K. Here true and test sin 2 26 13 = 0.1. NH is true and IH 
is test. Different curves correspond to various true 5c p in lower half plane (left panel) and upper half plane 
(right panel). 



From the right panel, we see that the degenerate (wrong hierarchy, wrong 5c p) solution for 
T2K occurs for 5 C p(test) around or ±180°. And in the range -120° < 5 C p (test) < -60°, 
where the degenerate wrong hierarchy solution for NOz/A occured, the x 2 f°r T2K is quite 
large. Because of this wide difference between the Sep values of the degenerate (wrong 
hierarchy, wrong Sep) solutions of NOz/A data and T2K data, together they rule out the 
wrong hierarchy. 

We illustrate this hierarchy discriminating power for a few cases where true value of 5c p is 
in the unfavourable half plane. Figs. 10, 11 and 12 show the x 2 plots for 5cp = 90°, 45° and 
respectively, with NH as the true hierarchy. The left panel shows x 2 f° r l-5*NOz/A alone 
whereas the right panel shows the x 2 f° r 1.5*NOz/A +2*T2K. These plots show x 2 f° r the 
two cases where the true and test hierarchies are the same and are opposite. In these plots, 
we have marginalised over sin 2 26 l 13 . In the left panel of Fig. 10, there is a large allowed 
region of 5cp(test) in the wrong half-plane, if the test hierarchy is the wrong hierarchy. In 
the right panel, this region is almost completely ruled out, with the addition of T2K data. 
There is a just a small region, centered around <5c*p(test) ~ 180°, where the x 2 dips just 
below 2.71, the cut-off for 90% C.L. We see very similar features for true 5cp = 45° in 
Fig. 11 and for true Sep = in Fig. 12. Essentially identical features are seen for the case 
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where IH is true hierarchy in Fig. 13 with true 6cp = —90°, Fig. 14 with true 5cp = —45° 
and Fig. 15 with true 5c p = 0. 



True (sin 2 2e 13 , 5 CP ) = (0.1, 90°) True (sin 2 29 13 , 8 CP ) = (0.1, 90°) 
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FIG. 10: (colour online) X 2 vs <5 C p(test) plots for 1.5*NOA (left panel) and 1.5*NOA + 2*T2K (right 
panel) with true sin 2 2#i3 = 0.1 and true Sep = 90°. 
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FIG. 11: (colour online) x 2 vs 5 C p(test) plots for 1.5*NO^A (left panel) and 1.5*NOA + 2*T2K (right 
panel) with true sin 2 26i% — 0.1 and true Sep — 45°. 
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True (sin 29 13 , 8 CP ) = (0.1,0) 



True (sin 26 13 , 8 CP ) = (0.1,0) 
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FIG. 12: (colour online) x 2 vs 5 CP (test) plots for 1.5*NO^A (left panel) and 1.5*NOi/A + 2*T2K (right 



panel) with true sin 2 2#i3 = 0.08 and true Sep — (systematics included). 



True (sin 2 29 13 , 5 CP ) = (0.1, -90°) True (sin 2 20 13 , 8 CP ) = (0.1, -90°) 
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FIG. 13: (colour online) X 2 vs <5 C p(tcst) plots for 1.5*NO^A (left panel) and 1.5*NO^A + 2*T2K (right 
panel) with true sin 2 2# 13 = 0.1 and true Sep — —90°. 
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True (sin 2 29 13 , 5 CP ) = (0.1, -45°) 



True (sin 2 20 13 , S CP ) = (0.1, -45°) 
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FIG. 14: (colour online) x 2 vs S C p(test) plots for 1.5*NOA (left panel) and 1.5*NOA + 2*T2K (right 




panel) with true sin 2 2#i3 = 0.1 and true Sep = —45°. 
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FIG. 15: (colour online) x 2 vs 8 C p (test) plots for 1.5*NOA (left panel) and 1.5*NOA + 2*T2K (right 

panel) with true sin 2 2(9i3 = 0.1 and true Sep = 0. 



Finally we consider how hierarchy sensitivity improves with increasing statistics. We 
consider three scenarios: 

• T2K will have a 5 year neutrino run with its design luminosity and NOz/A will run 
according to its present plan. 

18 



• T2K will have twice the above statistics and NOi^A will have 1.5 times its designed 
statistics. 



• T2K will have four times the above statistics and NOz/A will have thrice its designed 
statistics. 

The exclusion plots are given in Fig. 16. For all points to the right of the contours, the 
wrong hierarchy can be ruled out. 
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FIG. 16: (colour online) Hierarchy exclusion plots for combined data from NO^A and T2K with various 

boosts in statistics when NH is true (left panel) and when IH is true (right panel). 



In the left panel we assumed NH is the true hierarchy and in the right panel we assumed IH 
is the true hierarchy. We see that increasing the statistics from nominal values to 1.5*NOi/A 
_l_ 2*T2K dramatically improves the ability to rule out the wrong hierarchy, if <5c"p(true) is in 
the unfavourable half-plane. Further improvement occurs if the statistics are increased even 
more. In particular, if sin 2 2^13 = 0.1 [11, 12] the hierarchy can be essentially established 
at 90% C.L., for any true value of Sep, with 1.5 times the designed statistics of NOz/A and 
twice the designed statistics of T2K. This point was noted previously in Ref. [46]. 

It is evident now that an experiment that can exclude the wrong 6c p plane effectively can 
be of great help in determining hierarchy when run in conjunction with NOi^A. We saw that 
T2K, with a short baseline and smaller matter effects, has such properties. We now inquire 
whether having an experiment with a baseline shorter than T2K, such as C2F, which is 130 
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km long, can help. For such a short baseline, P^e is maximum at E = 0.25 GeV. At such 
energies, the matter term A is very small. 

To make a just comparison in terms of cost, we assume C2F to have the same beam 
power and detector size as that of T2K and 3 years each of v and v running. We consider 
two scenarios. NOz/A with 1.5 times its designed statistics and T2K with twice its designed 
statistics (scenario A) against NOz/A with 1.5 times its designed statistics and T2K and C2F 
with their nominal designed statistics (scenario B). In Fig. 17, we compare the ability of 
scenario A (left panel) and scenario B (right panel) to exclude the wrong hierarchy - wrong 
5cp region. The two panels are essentially identical. We found that scenarios A and B give 
the same allowed regions for all true values of Sep in the unfavourable half plane. Therefore, 
a shorter baseline experiment (L ~ 130 km) will not help in hierarchy determination. 



1.5*NOvA+2*T2K, 90% C.L. 



1.5*NOvA+T2K+C2F, 90% C.L. 



u 
to 



180 
120 
60 

-60 
-120 



-180 





, , , 




, , , 


' ' ' . 






X 










\J 










• i 

V 








Tru 
Tr 
. . . 


true point 

e:J>*H/Test:NH 
j4:NH«/ Test:IH 
«. . I i . . . 


x j 
. . . ■ 



0.05 0.07 0.09 0.11 0.13 0.15 

sin 2 20j 3 (test) 



0. 

u 



180 
120 
60 

-60 
-120 



-180 



' ' ' 


, , , 




, , , 


' ' ' . 






X 


























" . . . 


Tru 
Tn 
. . . 


true point 
e:NH / Test:NH 
je;N^ / TestrlH 
! . '.. i . . . 


x ■ 

. . . ■ 



0.05 0.07 0.09 



0.11 



0.13 0.15 



sin 2 29 1 3 (test) 



FIG. 17: (colour online) Allowed sin 2 26 13 -6 C p plots for 1.5*NOz/A + 2*T2K (left panel) and 1.5*NOi/A 
+ T2K + C2F (right panel) with true sin 2 28 13 = 0.1 and true S CP = 90°. 



IV. MEASURING 5 C p WITH P^ e 

A. 5qp measurement with T2K alone 

In the previous section, we discussed the capability of NOz/A and T2K to determine 
the mass hierarchy. We now turn our attention to the measurement of 8c p- Because of 
the hierarchy-5cp degeneracy, the determination of these two quantities go hand in hand. 
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Matter effects, which are hierarchy-dependent, induce a CP-like change in the oscillation 
probabilities. Therefore, it is expected that the effects of these two parameters can be 
disentangled by choosing baselines and energies where matter effects are small. Thus, a 
natural choice for accurate measurement of 5c p seems to be an experiment with a short 
baseline and low energy, like T2K or C2F. But, here we demonstrate that 5c p can not be 
measured in such experiments without first determining the hierarchy. For the purpose of 
this demonstration, in this subsection alone, we will assume that T2K will have equal three 
year runs in neutrino and anti-neutrino modes. This is done because such runs have the 
best capability to determine 5c p- However, even in such a case, 5c p can't be determined 
without first determining the hierarchy. 

In the following, we present 'allowed Sep' graphs. In generating these, we have kept 
sin 2 26'i3 fixed at 0.1. The graphs are plotted in the true 5cp-test 5cp plane. For every 
true value of 5cp, we indicate the range in test 5cp that can be excluded at 90% C.L. The 
plots have been shown for both true and wrong hierarchies. The dotted range, defined by 
X 2 < 2.71, shows the values of test 5cp that are compatible with the data, generated with 
5c p (true) as input. For a given true value of 5cp, the error in measuring 5c p is indicated 
by the spread of the dotted range along that 5cp{tme) vertical line. 



True : NH/Te st : NH , 90% C.L. True:NH/Test:IH, 90% C.L. 




hierarchy is normal (left panel) and inverted (right panel). 

Figure (18) shows the allowed 5cp plot for T2K. The points on the thick dashed line 
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in this figure correspond to the values of (fop (test) for which x 2 is minimum. If the test 
hierarchy is the same as the true hierarchy, then the \ 2 minimum occurs for Sep (test) = 
<5cp (true) and the allowed range of test Sep is around true Sep- But, if the test hierarchy 
is the wrong hierarchy, then the minimum of x 2 occurs for Scp{test) 7^ <5cp(true) and, in 
general, these two points are widely separated. This already gives a hint that an accurate 
measurement of Sep is not possible without first determining the hierarchy. This point is 
made more dramatic, when we consider the situation with more data from T2K. Fig. 19 
shows the allowed Sep plot for 10 times the statistics of T2K. For Sep (true) in the middle 
of the favourable half plane (—140° < Sep < —40°), the wrong hierarchy solution is ruled 
out. Thus both the hierarchy and the correct range of Sep are simultaneously determined. 
For all other values of Sep (true), we get a wrong value of Sep, if we assume the wrong 
hierarchy. For example, we see from the right panel of Fig. 19, for true Sep = —30°, we find 
that —130° < <5(7p(test) < —70°, when the test hierarchy is the wrong hierarchy. Similarly 
for true S CP = +60°, we find 140° < 5 CP (test) < 200°(= -160°). In particular, if true S CP is 
—10°, close to the CP conserving value 0, we have —150° < <5<7p(test) < —50°, encompassing 
maximal CP violation. The situation is similar for true Sep = —170°. Conversely, for true 
Sep — 90°, we have two allowed regions between to 40° and 140° to 180°, both of which 
are close to CP conservation. This figure makes it clear that it is impossible to have a 
measurement of Sep if we do not know the correct hierarchy. In fact, we are likely to get a 
completely misleading estimate of Sep if we assume the wrong hierarchy. The corresponding 
figures for C2F experiment show similar features. 
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B. 5c p measurement with T2K and NOz^A 

In this subsection, we consider the Sep measuring capability of NOz/A and T2K together. 
Here we revert back to the original assumption that T2K will run in neutrino mode only 
for 5 years. Fig. 20, shows the allowed 5cp plot of NOz/A , assuming NH is true. If the test 
hierarchy is the true hierarchy, the allowed range of 5c p will surround true Sep- If the test 
hierarchy is the wrong hierarchy we obtain a large allowed range with Sep far from the true 
value. 
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True : NH/Test : NH , 90% C.L. True:NH/Test:IH, 90% C.L. 




Test hierarchy is normal (left panel) and inverted (right panel). 

Fig. 21 shows the allowed 8c p plot for NOi^A and T2K together. In the left panel, the 
allowed range 5c p for the true hierarchy is shown. We see that this range is mostly in the 
correct half-plane. For wrong hierarchy, shown in the right panel, the large allowed region 
in wrong half plane is reduced, but a substantial region is still allowed. For the case where 
IH is the true hierarchy, similar features occur. 
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True : NH/Te st : NH , 90% C.L. True:NH/Test:IH, 90% C.L. 




0.1. Test hierarchy is normal (left panel) and inverted (right panel). 

If the statistics are increased to 1.5*NOz/A + 2*T2K, as seen in Fig. 22 then most of 
the wrong hierarchy allowed region is ruled out as already noted in section 3. For the true 
hierarchy, the allowed region is centered around true 5c p and is mostly in the correct half- 
plane. For the CP conserving case 5c p = (5c p = ±180°), there is a small additional 
allowed region around 5c p = ±180° (5c p = 0) but for which \ 2 is higher. If we limit our 
attention to the regions around Xmini then 1.5*NOzM. + 2*T2K can measure 5cp with an 
accuracy of ±40° for true 5c p = and ±60° for true 5c p = ±90°. 

It is curious that the CP conserving values of 5cp can be measured with better accuracy 
than large CP violating values. However, this point can be understood very simply in terms 
of Eq. (3). 5cp occurs in this equation as cos(A + 5cp)- Any experiment is designed such 
that the flux peaks at the energy where A rs 90°. Thus the 5c p term is approximately 
— sin 5cp- The slope of sin a; is large at x ~ or 180° and is very small at x ~ ±90°. 
Therefore the uncertainty in 5c p is small near or 180° and is large when 5c p is close to 
±90°. 

Thus we are led to the following important conclusion: 1.5*NOvA + 2*T2K can essen- 
tially determine the hierarchy and also give an allowed region of 5c p centered around its 
true value. Doubling of statistics will not lead to too much improvement in the allowed 
range of 5c p- Further strategies are needed to measure 5c p to a good accuracy. 
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A recent paper [48] envisaged some future very long baseline superbeam experiments. 
They found that the early data from these will determine hierarchy, and additional data is 
needed to measure Sep- We find that in the current scenario also, these considerations hold 
true. 



True : NH/Test : NH , 90% C.L. True:NH/Test:IH, 90% C.L. 
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FIG. 22: (colour online) Allowed S C p plots for 1.5*NOA + 2*T2K. Here NH is true. True and test 
sin 2 2#i3 = 0.1. Test hierarchy is normal (left panel) and inverted (right panel). 



V. SUMMARY 



In this paper we explored the hierarchy - 5c p degeneracy of P^ e of medium long 
baseline experiments. This degeneracy severely limits the ability of any single experiment 
to determine these quantities. The observed moderately large value of 9i% is certainly 
a very good news for the upcoming NOi^A, as it will lie in the region where NOi^A has 
appreciable reach for hierarchy determination if the value of 5cp happens to be favourable. 
We define the concept of favourable half-plane of Sep and show that the LHP(UHP) is 
the favourable (unfavourable) half-plane for NH and vice- verse for IH. We also show that 
NOi^A by itself can determine the hierarchy if 5qp is in the favourable half-plane and 
sin 2 2^i3 > 0.12. When 5qp is in the unfavourable half-plane, the data from NOz/A and 
T2K beautifully complement each other to rule out the wrong hierarchy. We explore the 
underlying physics in detail and deduce the statistics needed for hierarchy determination. 
Given the current best fit of sin 2 2# 13 ~ 0.1, the combined data from NO^A and T2K can 
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essentially resolve mass hierarchy for the entire 5cp range if the statistics for NOz/A and 
T2K are boosted by factors 1.5 and 2 respectively. A baseline of ~ 130 km will not be a 
bonus, over and above T2K, unless supplemented by huge statistics. 

In the last section we estimate the 5c p reach of NOz/A and T2K. We demonstrate that 
without knowing the hierarchy, measuring 5cp would be impossible. With 1.5*NOz/A + 
2*T2K, the allowed region of 5 C p is centered around its true value and is mostly in the correct 
half-plane. Here also, a short baseline of ~ 130 km will not provide better information than 
T2K with the same statistics. 
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